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The molecular mechanisms underlying specification from embryonic stem cells (ESCs) and maintenance of neural 
progenitor cells (NPCs) are largely unknown. Recently, we reported that the Zuotin-related factor 1 (Zrfl) is 
necessary for chromatin displacement of the Poly comb-repressive complex 1 (PRC1). We found that Zrfl is 
required for NPC specification from ESCs and that it promotes the expression of NPC markers, including the key 
regulator Pax6. Moreover, Zrfl is essential to establish and maintain Wnt ligand expression levels, which are 
necessary for NPC self -renewal. Reactivation of proper Wnt signaling in Zrfl -depleted NPCs restores Pax6 
expression and the self -renewal capacity. ESC-derived NPCs in vitro resemble most of the characteristics of the 
self-renewing NPCs located in the developing embryonic cortex, which are termed radial glial cells (RGCs). 
Depletion of Zrfl in vivo impairs the expression of key self-renewal regulators and Wnt ligand genes in RGCs. 
Thus, we demonstrate that Zrfl plays an essential role in NPC generation and maintenance. 
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Epigenetic regulators, such as proteins of the Polycomb- 
repressive complex 1 (PRC1) and PRC2, repress develop- 
mental genes and play a crucial role in stem cells, in line 
with their essential role in embryonic development 
(Schuettengruber and Cavalli 2009; Aloia et al. 2013). 
We showed recently that the chromatin component 
Zuotin-related factor 1 (Zrfl) derepresses developmental 
genes by displacing the PRC1 complex from chromatin 
during differentiation of human teratocarcinoma cells 
(Richly et al. 2010). Furthermore, the Caenorhabditis 
elegans Zrfl ortholog dnjl 1 was implicated in asymmetric 
cell division of neurosecretory motoneuron neuroblasts 
(Hatzold and Conradt 2008). Although these studies sug- 
gest that Zrfl is possibly involved in embryonic develop- 
ment, whether it is involved in stem cell maintenance and 
differentiation has not yet been explored. 

Mouse embryonic stem cells (mESCs) derived from 
embryonic day 3.5 (E3.5) blastocysts are a valid tool to 
study embryonic development since they are able to give 
rise to all cell types of the embryo in vivo and in vitro 
(Boiani and Scholer 2005; Niwa 2007; Nichols and Smith 
2012). In the last few years, several experimental methods 
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to efficiently generate neural progenitor cells (NPCs) 
from ESCs have been described (Conti et al. 2005; 
Colombo et al. 2006; Bibel et al. 2007). The first NPCs 
in the embryo are neuroepithelial cells present in the 
neural tube at El 1.5. Neuroepithelial cells give rise to 
Pax6-positive radial glial cells (RGCs) that remain at the 
ventricular zone (VZ) of the embryonic cortex (Campbell 
and Gotz 2002; Malatesta et al. 2003, 2008; Gotz and 
Huttner 2005; Kriegstein and Alvarez-Buylla 2009). RGCs 
act as self-renewing cellular elements, undergoing sym- 
metric proliferative divisions to maintain the pool of 
progenitors. RGCs also undergo asymmetric neurogenic 
division to generate more differentiated progenitors and 
to support migration of differentiating progenitors by 
acting as cellular scaffolding units. In adults, NPCs that 
originated from embryonic RGCs are present in restricted 
areas of the central nervous system (CNS), where they 
produce differentiated neurons and glia (Temple 2001; 
Merkle et al. 2007; Zhao et al. 2008; Miller and Gauthier- 
Fisher 2009; Fuentealba et al. 2012). Importantly, NPCs 
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are capable of repairing brain injuries upon transplanta- 
tion (Temple 2001; Aboody et al. 2011). 

Several transcription factors, including Sox 1/2, Neurog 
1/2, and Brn 1/2, are involved in neuroectodermal differ- 
entiation and specification of NPCs (Lee 1997; Aubert 
et al. 2003; Bani-Yaghoub et al. 2006; Suter et al. 2009; 
Gomez-Lopez et al. 2011). Moreover, the transcription 
factor Pax6 has been reported to be crucial for neural 
development and embryonic and adult NPC features, 
such as proliferation and differentiation (Stoykova et al. 
1996; Gotz et al. 1998; Heins et al. 2002; Sansom et al. 
2009; Tuoc et al. 2009). Several signaling pathways, such 
as Notch, Sonic hedgehog, and Wnt, are involved in the 
specification and maintenance of NPCs (Mizutani and 
Saito 2005; Han et al. 2008; Nusse 2008; Alvarez-Medina 
et al. 2009; Bluske et al. 2012; Bowman et al. 2013). 
During canonical Wnt signaling, (B-catenin accumulates 
in the cytoplasm before translocating to the nucleus, 
where it can activate transcription (MacDonald et al. 
2009). Wnt activity has been reported to antagonize 
the generation of NPCs from ESCs (Aubert et al. 2002; 
Haegele et al. 2003; Cajanek et al. 2009; Blauwkamp et al. 
2012), although p-catenin was found to be required for 
neural differentiation (Otero et al. 2004; Lyashenko et al. 
2011). Importantly, Wnt signaling is required for NPC 
maintenance (Chenn and Walsh 2002; Kalani et al. 2008) 
and differentiation into neurons (Hirabayashi et al. 2004; 
Lie et al. 2005). Recent evidence has shown that endog- 
enously produced Wnt ligands are important for self- 
renewal and multipotency of NPCs (Wexler et al. 2009). 
However, molecular mechanisms regulating Wnt ligand 
expression in NPCs are largely unknown. 

Here we show that Zrfl is essential for neuroectoder- 
mal specification and that it is required for NPC self- 
renewal by regulating Pax6 expression. Moreover, we 
unveil a role for Zrfl in the control of Wnt ligand expression 
in NPCs. We also demonstrate that in vivo Zrfl depletion in 
the embryonic VZ of the cortex impairs the expression of 
key self -renewal regulators. We propose that Zrfl is a key 
player required for first inducing NPC specification from 
ESCs and then maintaining NPC identity. 

Results 

Zrfl is required for neuroectodermal specification 
of ESCs 

To investigate whether Zrfl has a role in ESCs, we stably 
depleted it from mESCs (Fig. 1 A). We found that it did not 
affect either cell growth or the expression of pluripotency 
markers, such as Oct4, Nanog, Sox2, Rexl, and alkaline 
phosphatase (Fig. IB; Supplemental Fig. S1A,B). To in- 
vestigate the differentiation potential of Zrfl -depleted 
ESCs, we injected shZrf 1 (Zrfl shRNA) or shCtrl (control 
shRNA) cells into SCID mice to induce teratoma forma- 
tion. Strikingly, immuno staining of teratomas derived 
from shZrfl ESCs revealed a significant reduction of 
neuroectodermal markers, such as Nestin, p3-tubulin, 
and GFAP (markers of NPCs, post-mitotic neurons, and 
astrocytes, respectively), whereas endodermal (FoxA2) 



and mesodermal (actin a smooth muscle [SMA]) markers 
were properly expressed (Fig. 1Q Supplemental Fig. S2A). 
In agreement with this, RT-qPCR analysis of teratomas 
indicated that Zrfl depletion impaired the expression of 
neuroectodermal markers, including Nestin, Pax6, and 
Soxl, whereas mesodermal and endodermal markers 
were expressed at the same level or slightly higher in 
Zrfl -depleted teratomas (Fig. ID). We thus decided to 
focus on the role of Zrfl in neuroectodermal lineage 
specification. 

To this end, we used an ESC line that expresses GFP 
under the control of the Soxl promoter (Soxl -GFP ESCs), 
which is one of the earliest neuroectodermal markers 
induced in the embryo (Suter et al. 2009). Soxl-GFP ESCs 
were cultured in serum-free medium supplemented with 
N2 and B27 (N2B27 medium), which is known to pro- 
mote neural induction (Ying et al. 2003). While control 
cells showed up-regulation of GFP expression, Zrfl de- 
pletion significantly decreased the percentage of GFP- 
positive cells at days 3 and 5 in N2B27 medium (Fig. IE). 
On the other hand, pluripotency markers were down- 
regulated in shZrfl cells cultured in N2B27 medium to 
the same level as in shCtrl cells (Supplemental Fig. S2B). 
Importantly, nonneural markers were not detectable in 
Zrfl -depleted cells grown in N2B27 medium at either 
the mRNA (data not shown) or protein (Supplemental 
Fig. S2C) level. In addition, Zrfl depletion impaired the 
expression of the early neuroectodermal marker Nestin 
(Fig. IF) and the marker of post-mitotic neurons, p3- 
tubulin (Fig. 1G). To demonstrate that the absence of Zrfl 
was responsible for the effects of the shRNAs, we in- 
troduced Flag- tagged human Zrfl cDNA into Soxl-GFP 
ESCs depleted for Zrfl, thus restoring correct expression 
of Zrfl. As expected, re-expression of Zrfl rescued the 
differentiation defect of shZrfl ESCs (Fig. 1E-G). 

To determine whether Zrfl is also involved in other 
differentiation processes, we cultured shZrfl cells in the 
presence of activin A for 5 d to induce the expression of 
mesodermal and endodermal markers. RT-qPCR analysis 
revealed that shZrfl ESCs did not block the up-regulation 
of several mesodermal and endodermal markers in this 
culture condition (Supplemental Fig. S2D), suggesting 
that Zrfl is specifically required for the ESCs to commit 
to neuroectodermal lineages. 

Zrfl is essential for the generation of RGC-like 
NPCs from ESCs 

We next analyzed whether Zrfl is involved in generating 
NPCs from ESCs by depleting Zrfl from Soxl-GFP cells 
that were then treated with retinoic acid (RA) to induce 
their differentiation into Pax6-positive, RGC-like NPCs 
(Fig. 2A ; Bibel et al. 2004, 2007). We found that depleting 
Zrfl significantly down-regulated Soxl (Supplemental 
Fig. S3A,B) but did not impair the expression of pluripo- 
tency markers as compared with control cells (Supplemen- 
tal Fig. S3C). To elucidate the molecular function of Zrfl 
during the generation of RGC-like NPCs, we performed 
gene expression profiling at day 6 of neural induction (see 
Supplemental Table 1). At this time point, Zrfl depletion 



GENES & DEVELOPMENT 183 




Ectoderm Endoderm Mesoderm 



Ectoderm 



Endoderm 



Mesoderm 



Serum free medium supplemented with N2 and B27 



shCtrl 



shZrfl 





— do 




dO 


dO 




d5 




d5 




d5 




i 1 

\ 








r\ 








1 





FITC(Soxl-GFP) 




day 5 



shCtrl 



shZrfl 



rescued 







_ 


68,0+8,0% 


9,4±5,3% 


66,3+11,4% 



day 5 



(53-tubulin 
Zrfl 
Flag 
Gapdh 



Figure 1. Zrfl is required for neuroectodermal specification of ESCs. [A, left) Western blot for Zrfl after infection of E14Tg2A ESCs with two 
different lentiviral shRNAs. {Right) Band intensity was quantified by using ImageJ software. Data represent the average of three independent 
experiments. (**) P < 0.01. (B) The relative mRNA expression of pluripotency markers after Zrfl depletion. Values are normalized to Gapdh. 
Data represent the average of three independent experiments. (C) Mean fluorescent intensity of immunostained cells in shCtrl or shZrfl 
teratomas quantified by using ImageJ software is depicted in the graph. Data represent the average of five selected teratomas (10 fields for each 
teratoma). [*)P< 0.05. (D) The relative mRNA expression of selected markers of three germ layers in entire teratomas. Values are normalized 
to Gapdh. The analysis was performed on five different teratomas. (*) P < 0.05. (E, left) FACS profile of Soxl-GFP expression ESCs grown in 
self -renewing condition (dO) and in N2B27 medium at day 5 (d5). {Right) The percentage of Soxl-GFP-positive cells at days 3 and 5 is reported 
in the graph. Data represent the average of three independent experiments. (**) P < 0.01. [F) Nestin immunostaining at day 5 in N2B27 
medium. Bar, 200 |xm. The percentage of positive (as compared with total DAPI [40,60-diamidino-2-phenylmdole]-positive) cells is reported in 
the picture. Data represent three independent experiments (seven fields were counted for each experiment). (G) Western blot for f33-tubulin in 
shCtrl, shZrfl, and rescued cells (coexpressing shZrfl and human Zrfl cDNA tagged with 3xFlag) at day 5 in N2B27 medium. 
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Figure 2. Zrfl is required for generation of RGC-like NPCs from ESCs. (A) Scheme of generation of RGC-like NPCs. Briefly, ESCs were 
grown as cell aggregates (CAs) in medium supplemented with 10% serum. At day 4, RA was added to the medium to promote the 
generation of Pax6-positive RGC-like NPCs. At day 8, CAs were disaggregated, and RGC-like NPCs were harvested for the analysis. [B] 
Heat map illustrating the expression changes of selected NPC markers in shCtri and shZrfl CAs at day 6 of neural induction. Each lane 
corresponds to an independent biological sample. (C) The relative mRNA expression of NPC markers during generation of RGC-like 
NPCs. Values are normalized to Gapdh. Data represent the average of three independent experiments. (*) P < 0.05; (* *) P < 0.01. (D) Zrfl 
ChlP-qPCR [left panel) and RingIB ChlP-qPCR [right panel) of selected NPC markers in ESCs (dO) and CAs at day 6 of neural 
specification (d6). Values are expressed as percentage of input. Data represent three different experiments. (**) P < 0.01. [E) Pax6 
expression at mRNA [left panel) and protein [right panel) levels during NPC specification from ESCs. The graph represents data from 
three different experiments. (*) P < 0.05; (**) P < 0.01. 
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up-regulated 1315 genes and down-regulated 1118 genes 
(Supplemental Fig. S3D). Gene ontology (GO) analysis 
revealed that the down-regulated genes were enriched 
for neuroectoderm markers and were mainly involved in 
transcription and embryonic development. On the other 
hand, a significant enrichment of markers of nonneural 
lineages was observed in the up-regulated gene set 
(Supplemental Fig. S3D). Specifically, Zrfl knockdown 
impaired the expression of several NPC key regulators, 
including Soxl, Neurogl, Newog2, and Neurodl (Fig. 2B). 
Moreover, Zrfl depletion affected the expression of 
markers of the forebrain (FoxGl), midbrain (including 
LmxlA and Pitx3), and hindbrain (including HoxB4 and 
HoxB8) (Supplemental Fig. S3E). As expected, Zrfl re- 
expression restored the expression of the mentioned genes 
(Fig. 2C). 

Chromatin-bound Zrfl increases with RA treatment 
(Richly et al. 2010), leading us to hypothesize that Zrfl is 
recruited to NPC genes upon ESC differentiation. Indeed, 
chromatin immunoprecipitation (ChIP) analysis revealed 
that Zrfl was recruited to the promoter regions of several 
NPC markers, such as Neurodl, Brnl, and Brn2 (Fig. 2D, 
left panel). To investigate whether Zrfl is involved in 
derepressing neural genes, we studied the occupancy of 
the PRC1 subunit RinglB at Zrfl-bound genes. ChIP 
experiments followed by RT-qPCR indicated that RinglB 
was bound to Zrfl target genes in ESCs, whereas RinglB 
occupancy decreased in differentiating cells (Fig. 2D, right 
panel). However, in Zrfl -depleted ESCs, RinglB occu- 
pancy at Zrfl targets was maintained also in differentiat- 
ing cells (Fig. 2D, right panel). Taken together, these data 
suggest that Zrfl is essential for PRC1 displacement on 
selected NPC markers and their transcriptional activa- 
tion upon differentiation. 

Finally, we found that Pax6 expression in shZrfl cells 
was slightly reduced at day 6 and further down-regulated 
at day 8 of NPC specification at both the mRNA and 
protein levels as compared with shCtrl cells (Fig. 2E). This 
indicated that Zrfl is necessary to generate Pax6-positive 
NPCs from ESCs. 

Zrfl regulates Pax6 expression and is essential 
for self-renewal of NPCs 

Since Zrfl was crucial for the up-regulation of several 
NPC markers during NPC specification and thus for their 
expression in NPCs, we hypothesized that it could also 
have a role in the maintenance of NPC features. To test 
this hypothesis, we generated self-renewing NPCs from 
wild-type ESCs that could be maintained for >50 passages 
in culture (Colombo et al. 2006). These self -renewing 
NPCs were then infected with Ctrl or Zrfl shRNAs. Zrf 1- 
depleted NPCs grew slower than shCtrl NPCs (Fig. 3A), 
and Zrfl depletion significantly decreased cell prolifera- 
tion as measured by BrdU incorporation (Fig. 3B). On the 
other hand, no changes in the percentage of apoptotic 
cells were observed after Zrfl knockdown (Supplemental 
Fig. S4A). At the molecular level, Zrfl depletion resulted 
in the down-regulation of several genes involved in 
maintaining NPC identify, such as Pax6, Emx2, Hes5, 



Nkx2.1, and Dlxl (Fig. 3C), but did not affect other NPC 
markers, such as Sox2 and SSEA1 (Fig. 3C ; Supplemental 
Fig. S4B). The specificity of the shRNAs in NPCs was 
verified by using a vector encoding the human Zrfl cDNA. 
Re-expression of Zrfl restored the correct levels of NPC 
markers in shZrfl NPCs (Supplemental Fig. S4C). 

The self-renewal capability of NPCs is characterized by 
sphere formation. Strikingly, Zrfl depletion significantly 
affected primary and secondary sphere formation (Fig. 3D; 
Supplemental S4D). Importantly, re-expression of cor- 
rect Zrfl levels restored the expression of NPC markers 
and overcame the self -renewal defect (Supplemental 
Fig. S4D), indicating that Zrfl is essential for self- 
renewal of NPCs. 

Since Pax6 is the key self-renewal regulator of embry- 
onic and adult NPCs (Sansom et al. 2009), we checked 
whether re-expression of Pax6 in shZrfl NPCs was able to 
restore proper self-renewal capacity. We expressed a Flag- 
tagged mouse Pax6 cDNA in shCtrl and shZrfl NPCs. 
Pax6 levels in shCtrl NPCs were slightly over-expressed 
with respect to the normal endogenous levels. Impor- 
tantly, overexpression of Pax6 in shZrfl NPCs restored its 
levels similar to that observed endogenously (Fig. 3E). 
Pax6-shZrfl NPCs properly formed primary and second- 
ary spheres, thus exhibiting correct self-renewal features 
(Fig. 3F,G). Thus, we conclude that Pax6 down-regulation 
was responsible for the self -renewal defect in Zrfl- 
depleted cells. Since Pax6 has been reported to be crucial 
for not only self-renewal but also neuronal differentiation 
(Gotz et al. 1998; Lie et al. 2005), we checked whether 
Zrfl -depleted NPCs could still differentiate toward post- 
mitotic neurons. We found that Zrfl depletion impaired 
the generation of p3-tubulin-positive neurons when NPCs 
were grown either in the absence of EGF and FGF2 
(Supplemental Fig. S4E) or with a B27 supplement (Sup- 
plemental Fig. S4F), indicating that Zrfl also plays a role 
in NPC differentiation. Pax6 re-expression significantly 
improved the differentiation capacity of shZrfl NPCs 
(Supplemental Fig. S4F). Taken together, these results 
indicate that Zrfl is required for both self-renewal and 
differentiation in NPCs. 

A PRC 1- dependent Zrfl involvement in neurogenesis 

To gain insight into Zrfl function, we next performed 
ChIP experiments followed by massive parallel sequenc- 
ing (ChlP-seq) of Zrfl in the two stem cell populations of 
ESCs and NPCs. We did not find significant Zrfl -binding 
sites in ESCs (—150 peaks, false discovery rate [FDR] = 
100%), whereas we found —12,000 significant peaks 
(FDR < 5%) corresponding to 5989 genes in NPCs (see 
Supplemental Table 2). GO analysis of Zrfl target genes 
showed that they were significantly enriched for neural 
lineages (Fig. 4A). To elucidate a possible PRC 1 -dependent 
role of Zrfl in neurogenesis at a genome-wide level, we 
took advantage of previously published RinglB ChlP-seq 
in ESCs (Morey et al. 2012) and performed RinglB ChlP- 
seq in NPCs. RinglB occupied 5535 genes in ESCs and 
>10,000 genes in NPCs (see Supplemental Table 2). As 
expected, RinglB target genes were enriched for different 
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Figure 3. Zrfl is required for self -renewal of ESC-derived NPCs. (A) Growth curve of shCtrl and shZrfl NPCs and rescued NPCs. Data 
are representative of three independent experiments. (*) P < 0.05; (**) P < 0.01. [B) FACS plot of BrDU staining after incubation of 20 
min in shCtrl and shZrfl NPCs. Data are representative of three independent experiments. (C) The relative mRNA expression of 
several NPC markers in shCtrl and shZrfl NPCs. Values are normalized to Gapdh. Data represent the average of three independent 
experiments. (*) P < 0.05; (**) P < 0.01. (D) Representative pictures of shCtrl and shZrfl NPCs grown as spheres. Bars, 400 |xm. (E, top 
panel) Western blot of Flag and Zrfl expression in shCtrl and shZrfl NPCs electroporated with Flag-tagged Pax6 cDNA. [Bottom panel) 
The relative mRNA expression levels of Pax6 in shCtrl and shZrfl after Pax6 overexpression. [F) Representative pictures of shCtrl and 
shZrfl NPCs after Pax6 overexpression grown as spheres. Bars, 250 |xm. (G, top panel) Number of shCtrl and shZrfl NPCs overexpressing 
Pax6 that gave rise to spheres over 1000 cells. {Bottom panel) Percentage of cells derived by disaggregation of primary spheres able to give 
rise to secondary spheres over the total number of cells. Data are representative of five different experiments. (**) P < 0.01. 



lineages, including mesodermal, neural, and embryonic 
lineages (Fig. 4A,B). Since most of the Ring IB-binding 
sites were found at the transcriptional start site (TSS), 
we analyzed only peaks at TSSs for both RinglB and 
Zrfl. 



Comparison of Zrfl -binding sites in NPCs and RinglB in 
ESCs indicated that common targets were mostly enriched 
for neural lineages (Fig. 4 A), suggesting that Zrfl is mainly 
recruited to PRC1 targets related to neurogenesis during 
NPC specification. Moreover, a significant enrichment 
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for neural lineages was also observed in PRC1 targets not 
bound by Zrfl (Fig. 4A), indicating that PRC1 also exerts 
a Zrfl -independent function in neurogenesis. By compar- 
ing common target genes of Zrfl and RinglB in NPCs, we 
found that most of the Zrfl targets colocalize with RinglB 
and were significantly enriched for neural lineages (Fig. 4B). 
In line with this, we found that RinglB occupancy was 



lower in NPCs as compared with ESCs on proneural 
genes (such as Dlkl and Olig2) that were bound by Zrfl 
in NPCs (Fig. 4C ; Supplemental Fig. S5A). Notably, RinglB 
occupancy decreased also on the proneural gene FoxGl, 
which was not bound by Zrfl (as also shown in Fig. 2D), 
confirming a Zrfl -independent role for PRC1 in neuro- 
genesis (Fig. 4C). On the other hand, RinglB occupancy 
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was not affected on nonneural genes (such as Gata4 
and Gata6) that were not bound by Zrfl (Fig. 4C). 
Finally, we observed that a subset of Zrfl targets was 
not bound by PRC1 in either ESCs or NPCs and that 
these were enriched for nonneural lineages (Supplemen- 
tal Fig. S5B), suggesting a possible PRC 1 -independent 
function of Zrfl in nonneural lineages. Taken together, 
these data highlight a PRC 1 -dependent role of Zrfl in 
neurogenesis. 

Zrfl controls the expression of selected Wnt 
ligands in NPCs 

Among the Zrfl targets in NPCs, we found several genes 
involved in different signaling pathways, such as Pdgfra, 
Dlkl, Dclkl, and several Wnt ligands. ChlP-qPCR anal- 
ysis confirmed that Zrfl was bound to these genes in 
NPCs (Supplemental Fig. S5C-E). Since the autocrine/ 
paracrine Wnt signaling generated by endogenous Wnt 
ligands plays a crucial role in controlling NPC self- 
renewal (Kalani et al. 2008; Wexler et al. 2009), we 
hypothesized that Zrfl could regulate Wnt ligand genes. 
We thus analyzed the expression levels of Zrfl -bound 
Wnt ligands in ESCs, NPCs, and neurons. We found that 
several Zrfl -bound Wnt ligands, such as Wntl, Wnt3, 
Wnt5b, and Wnt6, were significantly up-regulated in 
NPCs as compared with ESCs. In addition, Wnt3 and 
Wnt6 were further up-regulated in neurons as compared 
with NPCs, whereas Wntlb was specifically up-regulated 
in neurons (Fig. 5A). Since Wnt signaling impairs NPC 
specification from ESCs ( Aubert et al. 2002; Haegele et al. 
2003; Cajanek et al. 2009; Blauwkamp et al. 2012), we 
analyzed the expression of active p-catenin at different 
developmental stages (e.g., in ESCs, during NPC specifi- 
cation, in self-renewing NPCs, and in post-mitotic neu- 
rons). In line with previous reports, the amount of active 
P-catenin decreased during NPC specification (Supple- 
mental Fig. S5F), indicating that the up-regulation of 
selected Wnt ligands in NPCs as compared with ESCs is 
not functional during NPC specification. To investigate 
whether Zrfl was directly involved in inducing the 
expression of the mentioned Wnt ligands, we performed 
Zrfl and RinglB ChIP in shCtrl and shZrfl NPCs. In- 
terestingly, Zrfl binding on all of the Wnt ligands 
analyzed significantly decreased in post-mitotic neurons 
as compared with NPCs, suggesting that Zrfl regulates 
Wnt ligand genes specifically in NPCs (Fig. 5B). As 
expected, we found that RinglB occupancy in control 
cells was significantly lower in NPCs and neurons as 
compared with ESCs. On the other hand, RinglB occu- 
pancy remained significantly higher in shZrf 1 NPCs (Fig. 
5C) as compared with control NPCs, suggesting that Zrfl 
plays an important role for RinglB displacement upon 
differentiation of ESCs into NPCs. Moreover, we found 
that several Zrfl -bound Wnt ligands were significantly 
down-regulated in Zrfl -depleted NPCs as compared with 
control cells at the mRNA and protein levels (Fig. 5D,E). 
Taken together, these data indicate that Zrfl is necessary 
to establish and maintain correct expression levels of 
selected Wnt ligands in NPCs. 



Regulation of Wnt signaling mediated by Zrfl controls 
NPC self-renewal 

To assess whether the down-regulation of Wnt ligands 
observed in Zrfl -depleted NPCs impaired the canonical 
Wnt signaling, we analyzed the amount of active p-catenin 
in shZrfl NPCs. Zrfl depletion resulted in decreased 
amounts of active p-catenin, indicating that canonical 
Wnt signaling was impaired in shZrfl NPCs (Fig. 6A). We 
therefore hypothesized that the impairment of the Wnt 
signaling observed in shZrfl NPCs could account for the 
self -renewal defects observed after Zrfl depletion. To test 
this hypothesis, shZrfl NPCs were treated with the 
GSK3 inhibitor CHIR99021 (CHIR), which is able to 
activate p-catenin in the absence of Wnt ligand stimula- 
tion (Fig. 6A ; Ying et al. 2008). Since our results indicated 
that down-regulation of Pax6 is responsible for the self- 
renewal defect of Zrfl -depleted NPCs (Fig. 3), we inves- 
tigated the level of Pax6 expression after CHIR treatment. 
Importantly, CHIR treatment at least partially rescued 
the expression level of Pax6 in shZrfl NPCs at both 
mRNA and protein levels (Fig. 6A,B), confirming that 
Pax6 is a downstream target of Zrfl in NPCs. In line with 
this, CHIR treatment resulted in a significant increase 
in the number of spheres generated by shZrfl NPCs 
(Fig. 6C). However, these spheres were smaller than 
control spheres, a finding compatible with the incom- 
plete rescue of Pax6 expression. CHIR treatment was 
able to completely rescue the capability of shZrfl NPCs 
to give rise to secondary spheres, thus restoring their self- 
renewal capacity (Fig. 6D). Partial rescue of Pax6 expres- 
sion and correct self -renewal features were also observed 
in shZrfl NPCs cultured with recombinant Wnt3a (Fig. 
6B-D) or R-spondin (data not shown). On the other hand, 
CHIR treatment was not able to restore the differentia- 
tion defect of shZrfl NPCs (Supplemental Fig. S5G). 
Taken together, these results unveil a novel role for 
Zrfl in Wnt signaling and indicate that it is essential 
for NPC self-renewal. 

Since Wnt signaling controls proliferation and dif- 
ferentiation of several stem cell types, we investigated 
a possible role for Zrfl in nonneural stem cell lineages. 
We performed Zrfl knockdown in primary mesenchymal 
stem cells (MSCs) (Supplemental Fig. S6A-C), the myo- 
blast cell line C2C12 (Supplemental Fig. S6D-F), and the 
ESC-derived hematopoietic progenitor cell line HPC7 
(Supplemental Fig. S7); these stem cell lineages are de- 
pendent on Wnt signaling. We observed that Zrfl de- 
pletion slightly decreased the growth rate of these stem 
cells, yet differentiation capability was not affected, 
suggesting that Zrfl exerts a specific role in NPCs. 

Zrfl regulates the expression of Wnt ligands 
and NPC master genes in vivo 

To investigate Zrfl function in NPCs in vivo, we in- 
troduced either the shCtrl or shZrfl vector with a GFP- 
expressing vector to cells in the VZ of an E13.5 cortex 
(Fig. 7A). At this developmental stage, the VZ contains 
Pax6-positive NPCs with RGC features. We then ana- 
lyzed the GFP-positive regions from El 5. 5 embryos (e.g., 
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Figure 5. Zrfl controls the expression of several Wnt ligands in NPCs. [A] The relative mRNA expression of Zrfl -bound Wnt ligands 
in ESCs, NPCs, and post-mitotic neurons. Data are representative of three different experiments. (**) P < 0.01. [B) Zrfl ChlP-qPCR of 
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48 h later). Interestingly, the number of shZrfl cells was 
significantly reduced as compared with shCtrl cells (Fig. 
7B ; Supplemental Fig. S8A). To test whether the reduc- 
tion of shZrfl cells could be due to a proliferation defect, 
we stained the cells with a Ki67 antibody, which marks 
actively proliferating cells. We found that Zrfl knock- 
down cells had fewer Ki67-positive cells (Supplemental 
Fig. S8B). On the other hand, we did not observe differ- 
ences in cell death, as a TUNEL assay showed no changes 
in shCtrl and shZrfl embryonic cortices between GFP/ 
shCtrl and GFP/shZrfl cells (Supplemental Fig. S8C). 



Finally, we investigated the migration capability of elec- 
troporated cells as an indication of the differentiation 
process. The number of GFP-positive cells able to migrate 
toward the proper cortical plate was severely reduced 
after Zrfl depletion (Fig. 7C), suggesting that Zrfl is 
necessary in RGC progeny to sustain a correct pattern of 
radial migration during cerebral cortex development. 

To gain insight into the Zrfl molecular function in the 
embryonic cortex, we sorted GFP-positive cells 48 h after 
electroporation with shCtrl- and shZrfl -encoding vec- 
tors. By analyzing these cell fractions by RT-qPCR, we 
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found that several markers of RGC identity, including 
Pax6 and Sox2, were down-regulated after Zrfl depletion 
at both the transcriptional and protein levels (Fig. 7D ; 
Supplemental Fig. S8D). Interestingly, Sox2 and Pax6 
expression was also diminished in surrounding nonelec- 
troporated cells (Supplemental Fig. S8D), suggesting that 
Zrfl depletion also exerts a non-cell-autonomous effect. 
Furthermore, in line with our data from cultured NPCs, 
mRNA levels of several Wnt ligands were down-regulated 
in shZrfl cells (Fig. 7D), suggesting that impairment of 
autocrine/paracrine Wnt signaling could be responsible 



for the non-cell-autonomous effect due to Zrfl depletion. 
Re-expression of Zrfl using a human Zrfl cDNA con- 
struct restored proper levels of Pax6 as well as of the other 
RGC markers, including Wnt ligand genes (Fig. 7D ; Sup- 
plemental Fig. S8D). Finally, we also analyzed the expres- 
sion of p3 -tubulin, a marker of post-mitotic neurons, in 
the electroporated regions of the embryonic cortex. Sur- 
prisingly, immunostaining analysis indicated that (33 -tu- 
bulin-positive neurons aberrantly accumulated in the VZ 
in Zrfl -depleted cells and their surrounding cells (Supple- 
mental Fig. S8E), suggesting that Zrfl depletion induced 
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Figure 7. Zrf 1 regulates the expression of Wnt ligands and self-renewal regulators in vivo. (A) Scheme of the experimental procedure. 
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premature differentiation in a significant fraction of 
RGCs. Taken together, these data highlight a role of Zrfl 
as key regulator of NPC features in vivo. 

Discussion 

Our data show a specific role of the chromatin-binding 
factor Zrfl in neural lineage specification and the self- 
renewal capability of Pax6-positive NPCs (summarized in 
Fig. 7E). Our earlier work had shown that Zrfl displaces 
the Polycomb proteins of the PRC1 complex from pro- 
moters upon differentiation of human teratocarcinoma 
lines (Richly et al. 2010). Here we show that Zrfl is spe- 
cifically involved in neural fate, as Zrfl depletion did not 
affect mesodermal and endodermal specification from 
ESCs or the stem cell features of several nonneural stem 
cell lineages. Furthermore, ChlP-seq data revealed that 
a PRC 1 -dependent role of Zrfl is required for neuro- 
genesis at a genome-wide level. Our data thus expand 
the molecular characterization of epigenetic regulation of 
neural development mediated by Polycomb proteins 
(Mohn et al. 2008; Hirabayashi et al. 2009; Aloia et al. 
2013). At the biochemical level, our previous data in 
human teratocarcinoma cells indicate that, in differentia- 
tion conditions, Zrfl can displace PRC1 from chromatin 
by competing for the binding to monoubiquitinated H2A. 
However, we cannot exclude that Zrfl recruitment to 
chromatin upon NPC specification actually promotes 
NPC generation, which eventually leads to the displace- 
ment of the PRC1 complex from its ESC targets. Impor- 
tantly, our ChlP-seq data suggest that Zrfl also exerts a 
Polycomb-independent role that should be explored more 
in detail in the future. 

At a molecular level, we show here that Zrfl sustains 
the expression of the key self-renewal regulator Pax6 and 
that re-expression of Pax6 restores the self-renewal capa- 
bility of Zrfl -depleted NPCs. Although several reports 
have elucidated the downstream targets of Pax6 (Sansom 
et al. 2009; Ninkovic et al. 2013; Walcher et al. 2013), so 
far, little is known about the regulation of Pax6 ex- 
pression. A possible mechanism is via the Wnt signaling 
pathway, which has also been reported to be crucial for 
NPC self-renewal (Chenn and Walsh 2002; Kalani et al. 
2008; Wexler et al. 2009). We found that Zrfl depletion 
induced a significant decrease in the amount of active 
P-catenin, the downstream effector of the canonical 
Wnt pathway. Treatment with the GSK3 -kinase inhibitor 
CHIR partially restored Pax6 expression and the self- 
renewal capacity of Zrfl -depleted NPCs, suggesting that 
Pax6 could be a downstream target of the Wnt signaling 
in NPCs. How does Zrfl depletion induce a decrease in 
the levels of active p-catenin? We found that Zrfl is 
essential for the establishment and maintenance of cor- 
rect levels of selected Wnt ligands in NPCs. Although 
autocrine/paracrine Wnt signaling sustained by endoge- 
nously produced Wnt ligands has been previously 
reported to be crucial for NPC self -renewal (Chenn and 
Walsh 2002; Kalani et al. 2008; Wexler et al. 2009), 
molecular mechanisms regulating the expression of 
Wnt ligands in NPCs have not yet been studied in detail. 



Our finding that Zrfl mediated the regulation of the 
Wnt ligands is in apparent contrast to the neurogenic 
function of Zrfl, as previous reports indicate that the Wnt 
signaling impairs neural induction from ESCs (Aubert 
et al. 2002; Haegele et al. 2003; Cajanek et al. 2009; 
Blauwkamp et al. 2012). However, in physiological con- 
ditions, we found that the levels of active p-catenin 
expressed were lower during neural specification than 
in either ESCs or self -renewing NPCs, suggesting that 
induction of selected Wnt ligands could be counteracted 
by specific inhibitors (Aubert et al. 2002; Haegele et al. 
2003; Cajanek et al. 2009). 

We also observed that Zrfl is required in vivo for self- 
renewal of embryonic neural progenitors (e.g., RGCs). 
Importantly, only few differences were observed in vitro 
and in vivo, the most evident of which was the significant 
down-regulation of Sox2 after Zrfl depletion in vivo. This 
discrepancy could be due to the high concentration of 
EGF in the NPC culture medium, which strongly sus- 
tains Sox2 expression and might compensate for any 
effects of Zrfl depletion in cultured NPCs (Hu et al. 
2010). The decreased level of Sox2 observed after intro- 
ducing shZrfl to the embryonic cortex could be respon- 
sible for the premature differentiation of RGCs into 
p3-tubulin-positive neurons, as previously described 
(Estivill-Torrus et al. 2002; Tiberi et al. 2012). On the 
other hand, in cultured NPCs, we observed a significant 
decrease in (33 -tubulin staining in conditions that pro- 
moted neuronal differentiation. We therefore do not be- 
lieve that p3 -tubulin up-regulation in shZrf 1 self -renewing 
RGCs is in contrast with our data regarding the proneuro- 
genic function of Zrfl in cultured NPCs. 

NPCs/neural stem cells have recently been used for 
tissue regeneration in model systems and are highly 
promising for developing treatments for neurodegenera- 
tive diseases (Temple 2001; Aboody et al. 2011). More- 
over, the recent finding that mouse embryonic fibroblasts 
(MEFs) can be reprogrammed into induced NPCs/neural 
stem cells has opened new horizons in regenerative 
medicine (Han et al. 2012; Lujan et al. 2012; Ring et al. 
2012; Thier et al. 2012). In this regard, our data suggest 
that, by developing strategies to recruit Zrfl to chromatin 
and sustain its transcriptional effect, Zrfl is a potential 
target for regenerative medicine. As Zrfl has been shown 
previously to be involved in oncogene-induced senes- 
cence (Ribeiro et al. 2013) and misregulated in several 
human tumors (Shoji et al. 1995; Resto et al. 2000; 
Greiner et al. 2003), we suggest that Zrfl could also play 
a role in brain tumor development. How Zrfl is recruited 
to chromatin and, once there, specifically regulates differ- 
ent processes (senescence, specification of NPCs, NPC 
maintenance, and differentiation) remains to be addressed. 
One possibility is that post-translational modifications 
of Zrfl could confer the specificity to Zrfl in different 
cell contexts, culture conditions, and developmental 
stages. Indeed, Zrfl is phosphorylated upon RA in P19 
cells (Smith et al. 2007). Another possibility is that the 
presence of specific interactors could specifically restrict 
Zrfl function at different developmental stages. For 
instance, the protein inhibitor of DNA binding (Idl) 
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has been shown previously to interact with Zrf 1 (Richly 
etal. 2010). 

Taken together, our data show that Zrfl is a crucial 
player involved in NPC specification. Initially, Zrfl acts 
through epigenetic processes to regulate the establish- 
ment of NPC fate. Once established, Zrfl continues to act 
via the Wnt signaling pathway and the neural transcrip- 
tional network to maintain NPC identity. 

Materials and methods 

Vectors 

The sequences of lentiviral shRNA #1/2 (plKO, Sigma) are listed 
in Supplemental Table S2. shCtrls were from Sigma. 

For rescue experiments in ESCs and NPCs, human Zrfl cDNA 
was cloned in the pCBA-3xFlag vector (Aloia et al. 2010). For 
rescue in the in vivo electroporation experiments, Zrfl was 
cloned in a pCAG-GFP vector (Sessa et al. 2008). For Pax6 
overexpression, mouse Pax6 cDNA (NM_0 13627) was cloned 
in pCBA-3xFlag vector. 

ESC culture 

E14Tg2A and Soxl-GFP ESC lines (a kind gift from Dr. Smith) 
(Ying et al. 2003) were cultured as previously described (Aloia 
etal. 2010). 

For alkaline phosphatase staining (Stemgent), 300 cells were 
plated per six-well plate and grown for 6 d. ESCs stably ex- 
pressing Zrfl shRNA #1/2/3 were generated by infection with 
lentiviral vectors and then selected with puromycin as pre- 
viously described (Morey et al. 2012). To obtain Zrfl -Flag ESCs, 
cells were transfected with pCBA Zrfl -3 X Flag vector using 
Lipofectamine 2000 according to the manufacturer's instruc- 
tions. Cells were then selected with neomycin for 7 d. 

ESC lineage specification 

For neuroectodermal specification, ESCs were cultured in se- 
rum-free medium supplemented with N2 and B27 as previously 
described (Ying et al. 2003). RGC-like NPCs were obtained by 
culturing ESCs as cell aggregates with RA as previously de- 
scribed (Bibel et al. 2004, 2007). Endodermal and mesodermal 
lineages were created by treating ESCs for 5 d with activin A (20 
ng/mL) as previously reported (Fei et al. 2010). Self -renewing 
NPCs were derived from ESCs as previously described (Colombo 
et al. 2006). 

NPC culture 

NPCs derived from ESCs were cultured in medium containing 
20 ng/mL EGF and 10 ng/mL FGF2 as previously described 
(Colombo et al. 2006). Zrfl and Pax6 re-expression was obtained 
by electropoating pCBA Zrfl 3xFlag and pCBA Pax6 3xFlag 
(Amaxa neural stem cell kit), and selected with neomycin for 
6 d. NPCs were treated with 3 |jlM CHIR (Deltaclon) or 100 ng/mL 
Wnt3a (R&D Systems) for 7 d, replacing fresh medium daily, as 
indicated. Sphere formation was obtained by plating 10,000 cells 
in a six-well plate (without adherent) for cell culture. Secondary 
spheres were generated by disaggregating primary spheres and 
plating a single cell in a 9 6- well plate. For differentiation toward 
post-mitotic neurons, NPCs were cultured in N2 medium with- 
out EGF or FGF2 for 5 d or with N2 medium supplemented with 
B27 without EGF or FGF2 for 7 d. 



MSC, HPC7, and C2C12 culture and differentiation 

HPC7 and C2C12 were cultured and differentiated as previously 
published (Pinto do O et al. 1998; Burattini et al. 2004). MSCs 
(Life SI 505 -100) were grown and differentiated as previously 
described (Meirelles Lda and Nardi 2003). 

ChIP assay 

ChIP assays were performed as previously described (Morey et al. 
2012) using anti-Zrfl and anti-RinglB antibodies (Richly et al. 
2010). Primers used for ChlP-qPCR analysis are listed in Supple- 
mental Table S3. Reads produced by ChlP-seq experiments of 
each of the samples were aligned with the mouse genome 
(version NCBIM37) using the Bowtie tool (version 0.12.7), 
allowing for two mismatches within the seed alignment. Se- 
quence tags aligned to the genome were subsequently analyzed 
with MACS software (version 1.4.1) to detect genomic regions 
enriched for multiple overlapping DNA fragments (peaks) that 
we considered to be putative binding sites. MACS estimated the 
FDR by comparing the peaks obtained from the samples with 
those from the control samples using the same P cutoff (1 X 10" 5 ). 
For this study, we considered all of the peaks <10% FDR for 
RinglB in ESCs, 1% FDR for RinglB in NPCs, and 5% FDR for 
Zrfl in NPCs. We used the Ensembl annotation (version 63) to 
annotate the peaks and detect the closest TSS. ChlP-seq signals 
around TSSs were calculated using BEDTools over TSS over- 
lapping peaks and normalizing the coverage for the total number 
of mapped reads. The Gene Expression Omnibus (GEO) acces- 
sion number for this set of data is GSE53542. 

Western blots 

Proteins were extracted as previously described (Aloia et al. 
2010). Antibodies used for Western blot were against Zrfl 
(Richly et al. 2010), Pax6 (1:250; Millipore), Gapdh (1:2000; Santa 
Cruz Biotechnology), (33-tubulin (1:300; Sigma), total (3-catenin 
(1:2000; BD Biosciences), active (3-catenin (1:1000; Millipore), 
Flag (1:500 ; Sigma), and Wntl (1:250; Sigma). 

RNA extraction and gene expression profile 

RNA was extracted with an RNA extraction kit (Qiagen). cDNA 
was generated from 1 mg of RNA with the First Strand cDNA 
synthesis kit (Fermentas). Primers used in the RT-qPCR assays 
are listed in Supplemental Table S3. For gene expression profiles, 
RNA from four independent experiments was isolated and 
processed as previously described (Morey et al. 2012). A cutoff 
of ± 2.0-fold with P < 0.01 was used to determine differentially 
expressed genes. The GEO accession number for this set of data 
is GSE53542. 

FACS analysis 

Cell suspensions were analyzed for GFP expression or stained 
with antibodies against SSEA1 (FAB2155A, R&D Systems), Mac-1 
(clone Ml/70, BD Biosciences), or Grl (clone RB6-8C5, BD Bio- 
sciences). Cells were analyzed with an LSR n FACS (BD Biosciences) 
using Diva version 6.1.2 (BD Biosciences) and Flowjo software 
version 10.0.6 

Cell immunofluorescence 

Cells were fixed in 4% paraformaldehyde and permeabilized 
with 0.2% Triton X-100 in 10% FBS (Invitrogen)/l% BSA in PBS 
for 15 min at room temperature. Samples were then incubated 
with primary antibodies overnight. The antibody used for cell 
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immunofluorescence was against (33-tubulin (1:300; Sigma) or 
eMHC (F1.652, Hybridoma Bank). Images were captured with an 
inverted microscope (DMI4000, Leica Microsystems). 

Teratoma formation 

SCID mice were injected subcutaneously with 5 X 10 5 ESCs 
infected with shCtrl or shZrfl #1 and #2. Four weeks after the 
injection, tumors were surgically dissected from the mice, 
homogenized (for total RNA was extracted) or embedded in 
OCT (Tissue-Tek) and frozen at — 80°C. Sections were stained 
with Nestin (1:300; Millipore), 03-tubulin (1:200; Covance), 
GFAP (1:300; Dako), SMA (1:200; Sigma), or Foxa2 (1:100; Abeam) 
antibodies. 

In vivo electroporation 

Electroporation was used to deliver expression vectors or re- 
porter constructs to the ventricular RGCs in utero as previously 
described (Saito and Nakatsuji 2001; Sessa et al. 2008, 2010). 
Briefly, uterine horns of E13.5 pregnant dams were exposed 
by midline laparotomy after anesthetization with 312 mg/kg 
Avertin. DNA plasmid ( 1 mL, corresponding to 3 mg) mixed with 
0.03% fast-green dye in PBS was injected into the telencephalic 
vesicle using a micropipette pulled through the uterine wall and 
amniotic sac. Platinum tweezer-style electrodes (7 mm) were 
placed outside the uterus over the telencephalon, and four pulses 
of 40 mV for 50 msec were applied at 950 msec intervals using 
a BTX square wave electroporator. The uterus was then placed 
back into the abdomen, the cavity was filled with warm sterile 
PBS, and the abdominal muscle and skin incisions were closed 
with silk sutures. 

Immunohistochemistry of embryonic cortices 

Tissues were fixed with 4% paraformaldeheyde overnight at 4°C, 
treated with 30% sucrose overnight at 4°C, and then set to the 
optimal cutting temperature for cryosectioning. Immunochem- 
istry analyses were performed as previously described (Colombo 
et al. 2004). Briefly, fixed cells or frozen sections were blocked in 
10% normal goat serum (Sigma) and 0.1% Triton X-100 (Sigma) 
for 1 h at room temperature. Cells or the sections were then left 
overnight at 4°C with the primary antibody in blocking solution. 
Secondary antibodies were applied to cells or sections for 2 h 
at room temperature. The following primary antibodies were 
used: chicken anti-GFP (1:500; Molecular Probes), rabbit anti-£3- 
tubulin (1:500; Covance), mouse anti-Sox2 (1:200; R&D Systems), 
and rabbit Pax6 (1:200; R&D Systems). Secondary antibodies 
(Molecular Probes) were conjugated with Alexa Fluor 488, 
Alexa Fluor 594, or Alexa Fluor 647 (1:500). DAPI (40,60- 
diamidino-2-phenylindole) was used for nuclear counterstaining. 
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